(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[      6587,        217]
NotebookOptionsPosition[      6053,        194]
NotebookOutlinePosition[      6390,        209]
CellTagsIndexPosition[      6347,        206]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"G", "[", 
    RowBox[{"r_", ",", "\[Theta]_"}], "]"}], " ", ":=", " ", 
   RowBox[{
    FractionBox["1", 
     RowBox[{"2", " ", "Pi", " ", 
      RowBox[{"\[Sqrt]", 
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["r", "2"], "+", "1", "-", 
         RowBox[{"2", "r", " ", 
          RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ")"}]}]}]], "+", 
    RowBox[{
     FractionBox["1", 
      RowBox[{"8", " ", "Pi", " "}]], 
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["r", "2"], "+", "1"}], ")"}]}], "+", 
    RowBox[{
     FractionBox["1", 
      RowBox[{"4", " ", "Pi"}]], 
     RowBox[{"Log", "[", 
      FractionBox["2", 
       RowBox[{"1", "-", 
        RowBox[{"r", " ", 
         RowBox[{"Cos", "[", "\[Theta]", "]"}]}], " ", "+", " ", 
        RowBox[{"\[Sqrt]", 
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["r", "2"], "+", "1", "-", 
           RowBox[{"2", "r", " ", 
            RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ")"}]}]}]], "]"}]}], 
    "-", 
    FractionBox["7", 
     RowBox[{"10", " ", "Pi"}]]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.515758767155343*^9, 3.515758839410973*^9}, {
  3.515758875707122*^9, 3.515758945664326*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Gr", "[", 
    RowBox[{"r_", ",", "\[Theta]_"}], "]"}], " ", ":=", 
   RowBox[{"D", "[", 
    RowBox[{
     RowBox[{"G", "[", 
      RowBox[{"r", ",", "\[Theta]"}], "]"}], ",", "r"}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.515758948639232*^9, 3.515758954596882*^9}, {
  3.51575900175814*^9, 3.515759017823431*^9}, {3.515760490194415*^9, 
  3.515760496658956*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Gr", "[", 
  RowBox[{"r", ",", "\[Theta]"}], "]"}]], "Input",
 CellChangeTimes->{{3.515758962894805*^9, 3.515758963840903*^9}, {
  3.515759024149865*^9, 3.515759031624698*^9}, {3.515759339518636*^9, 
  3.515759361785986*^9}, {3.515759802044888*^9, 3.515759806167813*^9}, {
  3.515760500747134*^9, 3.515760522336472*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["r", 
   RowBox[{"4", " ", "\[Pi]"}]], "-", 
  FractionBox[
   RowBox[{
    RowBox[{"2", " ", "r"}], "-", 
    RowBox[{"2", " ", 
     RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], 
   RowBox[{"4", " ", "\[Pi]", " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"1", "+", 
       SuperscriptBox["r", "2"], "-", 
       RowBox[{"2", " ", "r", " ", 
        RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ")"}], 
     RowBox[{"3", "/", "2"}]]}]], "-", 
  FractionBox[
   RowBox[{
    RowBox[{"-", 
     RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "+", 
    FractionBox[
     RowBox[{
      RowBox[{"2", " ", "r"}], "-", 
      RowBox[{"2", " ", 
       RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], 
     RowBox[{"2", " ", 
      SqrtBox[
       RowBox[{"1", "+", 
        SuperscriptBox["r", "2"], "-", 
        RowBox[{"2", " ", "r", " ", 
         RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}]]}]]}], 
   RowBox[{"4", " ", "\[Pi]", " ", 
    RowBox[{"(", 
     RowBox[{"1", "-", 
      RowBox[{"r", " ", 
       RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "+", 
      SqrtBox[
       RowBox[{"1", "+", 
        SuperscriptBox["r", "2"], "-", 
        RowBox[{"2", " ", "r", " ", 
         RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}]]}], ")"}]}]]}]], "Output",
 CellChangeTimes->{
  3.515759031954256*^9, {3.515759349977931*^9, 3.515759362189228*^9}, 
   3.515759806574009*^9, {3.515760502665165*^9, 3.515760523159166*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"K", "[", 
    RowBox[{"r_", ",", "\[Theta]_"}], "]"}], ":=", 
   RowBox[{
    FractionBox["r", 
     RowBox[{"4", " ", "\[Pi]"}]], "-", 
    FractionBox[
     RowBox[{
      RowBox[{"2", " ", "r"}], "-", 
      RowBox[{"2", " ", 
       RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], 
     RowBox[{"4", " ", "\[Pi]", " ", 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"1", "+", 
         SuperscriptBox["r", "2"], "-", 
         RowBox[{"2", " ", "r", " ", 
          RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ")"}], 
       RowBox[{"3", "/", "2"}]]}]], "-", 
    FractionBox[
     RowBox[{
      RowBox[{"-", 
       RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "+", 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "r"}], "-", 
        RowBox[{"2", " ", 
         RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], 
       RowBox[{"2", " ", 
        SqrtBox[
         RowBox[{"1", "+", 
          SuperscriptBox["r", "2"], "-", 
          RowBox[{"2", " ", "r", " ", 
           RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}]]}]]}], 
     RowBox[{"4", " ", "\[Pi]", " ", 
      RowBox[{"(", 
       RowBox[{"1", "-", 
        RowBox[{"r", " ", 
         RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "+", 
        SqrtBox[
         RowBox[{"1", "+", 
          SuperscriptBox["r", "2"], "-", 
          RowBox[{"2", " ", "r", " ", 
           RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}]]}], ")"}]}]]}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.515759034267208*^9, 3.515759048175124*^9}, {
  3.515759168744759*^9, 3.515759172099561*^9}, {3.515759203659821*^9, 
  3.515759216365253*^9}, {3.515760593977547*^9, 3.515760606206832*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", 
  RowBox[{"K", "[", 
   RowBox[{"1", ",", "\[Theta]"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.515760610039792*^9, 3.515760618891164*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{3.51576061916284*^9}]
}, Open  ]]
},
WindowSize->{740, 607},
WindowMargins->{{118, Automatic}, {Automatic, 0}},
FrontEndVersion->"8.0 for Linux x86 (32-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 1250, 38, 111, "Input"],
Cell[1810, 60, 422, 12, 30, "Input"],
Cell[CellGroupData[{
Cell[2257, 76, 344, 6, 30, "Input"],
Cell[2604, 84, 1448, 44, 80, "Output"]
}, Open  ]],
Cell[4067, 131, 1683, 49, 152, "Input"],
Cell[CellGroupData[{
Cell[5775, 184, 190, 4, 30, "Input"],
Cell[5968, 190, 69, 1, 30, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)

